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A .  Problems  Studied 


The  major  problem  considered  in  this  grant  was  the 
investigation  oT  the  mathematical  properties  of  a  new  class  of 
oscillators.  These  oscillators  have  the  property  that  their 
defining  differential  equation  can  become  unbounded  for  finite 
values  of  the  dependent  variable  or  its  derivative.  We  have 
called  such  systems,  “singular,  nonlinear  oscillators.”^ 
Examples  of  such  systems  occur  in  electronic  circuits,^  plasma 
physics,  and  nonlinear  vibrations.  ’ 

A  model  differential  equation  that  has  the  required 
mathematical  properties  is  the  so-called  WCM  equation® 


A*  —  . 

X  +  X  =  £  7 - 5  X , 

1  -  x-* 


where  £  and  /i  a^re  positive  parameters.  Of  interest  is  the  fact 
that  the  usual  methods  (perturbation  techniques,  multi-time 
expansions,  method  of  slowly  varying  amplitude  and  phase),  that 
depend  on  e  being  small,  cannot  be  applied  to  determine 
analytic  approximations  to  the  solutions  of  Eq .  (1).^  Hence, 
our  investigations  have  centered  on  the  following  questions: 

(i)  What  are  the  general  solution  behaviors  of  Eq .  (1)? 
In  particular,  what  can  be  stated  about  the  possible  existence 
of  limit-cycles,  their  number,  uniqueness,  location,  etc.? 

(ii)  For  the  construction  of  approximate  analytic 
solutions,  what  should  be  the  expansion  parameter?  Its  clearly 
not  (I  Is  it  a  function  of  /i?  If  so,  what  function  of  /i? 

(iii)  Equation  (1)  can  be  considered  a  large  amplitude^’^ 
version  of  the  van  der  Pol  equation^ 


X  +  X  =  £(1  -  x*)x. 


See  Bibliography  for  references. 
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What  is  "the  relationship  between  their  solutions?  How  do  their 
phase-spaces  relate  to  each  other?  What  is  the  proper  limit 
for  Eq.  (1)  to  reduce  to  Eq.  (2)?  It  is  clear  that  the  “e”  in 
the  two  equations  are  not  the  same- 

(iv)  The  limit-cycle  of  the  van  der  Pol  equation  can  be 
investigated  using  the  Hopt  biTurcation  theorem.  Can  the  same 
be  done  for  Eq.  (1)? 

(v)  Simple  harmonic  balance  allows  the  calculation  of  an 
excellent  analytic  approximation  to  the  solution  of  Eq .  (1)-^ 
What  is  the  best  procedure  for  obtaining  higher-order 
approximations  using  generalized  harmonic  balance  techniques? 

(vi)  What  new  phenomena  occur  when  two  or  more  singular 
oscillators  are  coupled?  Do  special  types  of  periodic 
solutions  exist?  If  so,  how  can  approximate  analytic  solutions 
be  obtained  and  their  stability  determined? 

The  following  section  B  provides  a  summary  of  the  most 
important  results  on  the  above  problems.  Section  C  gives  a 
list  of  all  abstracts  and  papers  that  have  been  or  will  be 
published  related  to  the  research  objectives  of  this  grant. 
Finally,  section  D  is  a  listing  of  all  senior  personnel  and 
graduate  students  who  have  been  involved  in  various  aspects  of 
this  research  project. 
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B .  Summary  of  Results 

The  following  is  a  concise  summary  of  "the  research 
accomplished  under  "this  grant.  These  results  have  been 
published  or  accepted  for  publication  in  refereed  scientific 
journals . 

References  in  this  section  refer  to  abstracts  or  papers 
given  in  section  C  of  this  report.  For  exeunple,  A2  and  P3  are, 
respectively,  the  second  abstract  and  the  third  paper  listed  in 
section  C. 

1  .  Failure  of  the  Method  of  Slowly  Varying  Amplitude 
and  Phase  [P8] 

In  paper  P8,  Mickens  and  Ramadhani  examined  the  reasons 
why  the  method  of  slowly  varying  amplitude  and  phase  cannot  be 
applied  to  singular,  nonlinear  oscillator  equations.  The 
breakdown  occurs  because  the  nonlinear  term  does  not  have  a 
uniform  bound  on  its  magnitude  in  arbitrary  finite  regions  of 
its  associated  phase-space.  Consequently,  the  defining 
equation  for  the  amplitude  has  singularities  that  do  not 
correspond  to  the  actual  dynamics  of  the  system. 

2.  Harmonic  Balance  [PI] 

Early  in  this  project,  we  showed  that  the  method  of  simple 
harmonic  balance  provides  an  easy  and  direct  procedure  for 
calculating  the  fundamental  properties  of  our  model  singular, 
nonlinear  oscillator.  First,  relations  giving  the  period  and 
the  amplitude  of  the  unique,  stable  limit-cycle  were  found. 
Secondly,  this  calculation  indicated  that  the  small  parameter  c 
was  not  the  relevant  expansion  parameter  for  constructing  a 
perturbation  type  solution. 
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3 .  A  Perturbation  Procedure  [A1 ,  P2,  P3] 

Since  e  is  not  the  parameter  that  can  be  used  for  the 
basis  of  a  perturbation  expansion,  what  parameter  should  be 
used?  We  approached  this  problem  by  considering  the  van  der 
Pol  oscillator  differential  equation.  The  reason  for  doing  so 
is  a  consequence  of  the  fact  that  our  model  singular,  nonlinear 
equation  can  be  considered  a  large  aunplitude  form  of  the  van 
der  Pol  equation. 

The  applications  of  the  Hopf  bifurcation  theorem  to  the 
van  der  Pol  equation  [P2]  lead  us  to  the  conclusion  that  the 
expansion  parameter  was  /<,  or  more  correctly  the  square-root  of 
ft  and  its  various  powers.  This  procedure  was  then  applied  to 
our  singular  equation  and  a  uniformly  valid  perturbation  series 
solution  was  found  [P3]  .  The  properties  of  our  equation, 
namely,  that  it  only  contained  “odd”  number  of  expressions  in 
its  rationalized  form,  allowed  for  a  greatly  simplified  final 
expression  for  the  perturbation  series  solution. 

4 .  Phase-Space  Techniques  [P4] 

An  advantage  of  examining  the  geometric  properties  of  the 
solutions  of  differential  equations  in  phase-space  is  that 
often  important  information  can  be  obtained  on  possible 
solution  behavior  without  much  initial  input  or  work.  In 
reference  P4 ,  we  showed  that  just  based  upon  phase-space 
considerations,  our  model  singular,  nonlinear  equation  had  to 
have  at,  least  one  stable  limit-cycle.  In  more  detail,  for 
€  >  0  and,  fi  >  0  and  small,  an  odd  number  of  limit-cycles  could 
exist.  The  application  of  the  Hopf  bifurcation  theorem  then 
showed  that  in  fact  only  one,  stable  limit-cycle  existed. 
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5 .  Transient  Behavior  [P5 ,  P73 

The  procedures  of  harmonic  balance  and  perturbation 
discussed  above  allow  only  the  determination  of  accurate 
approximations  to  the  steady-state  of  our  model  singular 
oscillator.  This  means  that  the  ajnplitude  and  angular 
frequency  can  be  calculated,  but,  not  the  transient  behavior  of 
the  solution  for  initial  conditions  not  on  the  limit-cycle.  In 
paper  P5 ,  we  showed  that  such  behavior  could  be  obtained  using 
a  variant  of  the  method  of  averaging.  The  stability  of  the 
limit-cycle  was,  as  a  by-product,  also  determined. 

There  are  many  dynamic  systems  in  engineering,  chemistry, 
biology  and  mathematics  that  can  be  modeled  in  terms  of  two 
coupled  first-order  ordinary  differential  qeuations.  In  paper 
P7,  we  extended  the  methods  of  paper  P5  to  cover  this  more 
generalized  situation.  [To  date,  I  have  received  more  than  80 
requests  for  reprints  of  this  paper!^ 

6 .  Coupled  Singular  Oscillators  [A2 ,  A5 ,  P9 ,  Pll] 

Coupled  nonlinear  oscillations  can  have  solutions  that  are 
periodic,  quasi -per iodic  or  chaotic.  The  possible  solution 
behaviors  not  only  depend  on  the  specific  values  of  the  system 
parameters  and  initial  conditions,  but,  also  on  how  the 
oscillators  are  coupled  to  each  other. 

Our  investigations  are  presented  in  two  abstracts,  A2  and 
A5 ,  and  two  papers,  P9  and  Pll.  The  particular  case  that  we 
studied  was  for  two  identical  singular  oscillators  interacting 
by  means  of  both  linear  and  nonlinear  couplings.  Our 
mathematical  analysis  used  the  concept  of  “nonlinear  normal 
modes.”  For  the  systems  studied,  we  found  that  there  were 
always  periodic  solutions  that  satisfied  one  of  the  symmetry 
cond it  ions 

x(t)  =  y(t)  or  x(t)  =  -y(t)  . 
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These  solutions  correspond,  respectively,  to  symmetric  and 
ant i -symmetric  nonlinear  normal  modes.  For  the  case  of  a 
particular  nonlinear  coupling,  the  symmetric  solution  was 
stable  and  the  anti -symmetric  solution  unstable.  However,  for 
the  linear  coupling  situation,  both  types  of  solutions  were 
stable.  Approximate  analytic  expressions  were  obtained  for 
these  solutions.  In  addition,  all  of  these  results  were 
verified  by  numerical  integration  of  the  equations  of  motion. 

It  should  also  be  stated  that  for  two  linearly  coupled 
singular  oscillators,  the  particular  “final”  periodic  state, 
that  the  system  approached  as  t  — ►  oo,  was  a  function  of  the 
choice  of  the  values  of  the  initial  conditions.  In  other 
words,  there  are  basins  of  attraction  for  the  symmetric  and 
ant i- symmetric  periodic  states. 

Finally,  evidence  was  found  for  so-called  “transition” 
states,  i.e.,  states  for  which  (numerically)  the  long-term 
behavior  appears  to  be  chaotic.  The  initial  conditions  for 
such  states  appear  to  be  localized  in  phase-space. 

7 .  Miscellaneous  Problems  [P6 ,  PIO] 

A  number  of  important  physical  systems  can  be  modeled  by 
van  der  Pol  type  differential  equations,  i.e., 

X  +  g(x)  =  ef(x)x. 


where 


g(x)  =  -g(-x),  f(x)  =  f(-x),  c  >  0. 


In  paper  P6 ,  we  exajnined  certain  mathematical  properties  of  the 
periodic  solutions  to  this  class  of  equations  as  a  function  of 
the  parameter  f. 
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A  very  interesting  nonlinear  oscillator  is  the  anti¬ 
symmetric  quadratic  oscillator,  i.e., 

X  +  X  +  f|x|x  =  0. 

We  showed  that  all  solutions  of  this  differential  equation  are 
bounded  and  periodic.  This  follows  directly  from  the  fact  that 
the  system  corresponds  to  a  conservative  oscillator  having  the 
potential  energy 


V(x) 


|x|x’ 


Analytical  approximations  to  the  solutions  have  been  calculated 
using  both  the  method  of  harmonic  balance  and  the  Krylov- 
Bogol iubov-Mitropolsky  procedure . 
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